1 



DISCRETE AND DIFFUSE X-RAY EMISSION IN THE NUCLEUS AND DISK OF THE 

STARBURST SPIRAL GALAXY M83 

R. Soria 1 and K. Wu 1 

Milliard Space Science Laboratory, University College London, Holmbury St Mary, Surrey RH5 6NT, UK 



(N 
O 
O 
(N 



> 

o 

(N 
O 
(N 
O 

6 

CO 



X 



Abstract 



et al. 1997). Twenty-one point sources were found in the 



We have studied the face-on, barred spiral M83 (NGC 
5236) with Chandra. Eighty-one point sources are detected 
(above 3.5-ct) in the ACIS S3 image: 15 of them are within 
the inner 16" region (starburst nucleus, resolved for the 
first time with Chandra), and 23 within the inner 60" (in- 
cluding the bar). The luminosity distributions of the sources 
in the inner 60" region (nucleus and stellar bar) is a sin- 
gle power law, which we interpret as due to continuous, 
ongoing star formation. Outside this inner region, there is 
a smaller fraction of bright sources, which we interpret as 
evidence of an aging population from a past episode of star 
formation. About 50% of the total emission in the nuclear 
region is unresolved; of this, about 70% can be attributed 
to hot thermal plasma, and we interpret the rest as due 
to unresolved point sources (eg, faint X-ray binaries). The 
unresolved X-ray emission also shows differences between 
the nuclear region and the spiral arms. In the nuclear re- 
gion, the electron temperature of the thermal plasma is 
w 0.58 keV. In the spiral arms, the thermal component is 
at kT » 0.31 keV and a power-law component dominates 
at energies > 1 keV. The high abundance of C, Ne, Mg, 
Si and S with respect to Fe suggests that the interstellar 
medium is enriched and heated by core-collapse supernova 
explosions and winds from massive stars. 
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1. Introduction 

M83 (NGC 5236) is a grand-design, barred spiral galaxy 
(Hubble type SAB(s)c) with a starburst nucleus. Distance 
estimates are still very uncertain. A value of 3.7 Mpc was 
obtained by |de Vaucouleurs et al. (1991) . This places the 
galaxy in the Centaurus A group, whose members have a 
large spread in morphology and high velocities, indicating 
that the group is not virialiscd and tidal interactions and 
merging are frequent ( |de Vaucouleurs 1979 ; |C6te et al 
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M83 wa s observed in the X-ra y bands by Einstein in 
1979- 1981 (iTrinchicri et al. 1985|) , by ROSAT in 1992- 
1994 Jlmmler et al. 1999|) , and by ASCA in 1994 flOkada 



ROSAT/HKl image, but the starburst nuclear region was 
unresolved. 

M83 was observed by Chandra on 2000 April 29, with 
the ACIS-S3 chip at the focus. The data became available 
to the public in mid-2001. The total exposure time was 
50.978 ks; after screening out observational intervals cor- 
responding to background flares, we retained a good time 
interval of 49.497 ks. 

In this paper we present the luminosity distribution of 
the discrete source population and discuss the properties 
of the unresolved emission in the nuclear region and in the 
disk. For further details on the data analysis techniques, 
and for more extensive discussions on the properties of the 
individual sources, see Soria fe Wu (2002) . 



2. Global properties of the discrete sources 

A total of 81 point sources are detected in the S3 chip at 
a 3.5-cr level in the 0.3-8.0 keV band. The source list is 



given in Soria & Wu (2002). Comparing the position of 



the Chandra S3 sources with a VLT B image shows that 
the off-centre sources tend to associate with the optically 
bright regions (Figure 1). The sources have a large spread 
in the hardness of their X-ray emission. A "true-colour" 
X-ray image of the nuclear region is shown in Figure 2, 
bottom panel. 

Separating the sources inside and outside a circular 
region of radius 60" from the geometric centre of the X- 
ray emission reveals that the two groups have different 
luminosity distributions in the 0.3-8.0 keV band. (A linear 
separation of 60" corresponds to 1.1 kpc for a distance of 
3.7 Mpc, and is roughly half of the total length of the 
major galactic bar.) The cumulative log N(>S) - log S 
distribution (where S are the photon counts) of the sources 
inside this inner region can be described as a single power 
law, with a slope of —0.8. The log N(>S) - log S curve 
of the sources outside the circular inner region instead 
is neither a single nor a broken power law (Figure 3). It 
shows a kink at S w 250 cts; the slope of the curve above 
the kink is —1.3, while it is —0.6 at the faint end. If we 
assume a foreground absorbing column density «h = 4 x 
10 20 cur 2 ( |5chlegel et al. 1998[ ), a distance of 3.7 Mpc and 



a power-law spectrum with photon index T = 1.5 for all 
the sources, 100 counts (« 2.0 x 10~ 3 cts s~ x ) correspond 
to an unabsorbed source luminosity L x — 2.3 xlO 37 erg s _1 
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Figure 1. Spatial distribution of the discrete X-ray sources 
overplotted on a VLT B image (size: 6'.8 x 6'.8). 



in the 0.3-8.0 keV band. The kink in the log N(>S) - 
log S curve of the sources outside the 60" circle is therefore 
located at L x « 6 x 10 37 erg s" 1 (0.3-8.0 keV band). 



We estimate from the Deep Field South survey (|Giac- 



coni et |al. 200l|) that about 15% of the 81 sources are back- 
ground AGN; the expected number in the inner 60" circle 
is smaller than one. The kink in the log N(>S) - log S 
curve for the outer sources and the values of the slope at 
both ends are unaffected by the background subtraction. 

The natter slope of the log N(>S) - log S curve at 
the high-luminosity end for the population of sources in 
the inner 60" region implies a larger proportion of bright 
sources than in source population further away from the 
nucleus. The situation is different for example in the spi- 
ral galaxy M81, where most bright sources are found in 



the galactic disk instead of the nuclear region (Tcnnant 
et al. S|001|)- If tne flatness of the slope in the log N(>S) 



log S curve is a characteristic of ongoing star formation 



(Wu 2001), the difference in the spatial distribution of the 
brightest sources in M83 and M81 is simply a consequence 
of the fact that M83 has a starburst nucleus while star for- 
mation in galaxies such as M81 is presently more efficient 
in the disk. 



Figure 2. Top panel: A low-resolution (ACIS-S3 chip, with 
pixels binned 4x4,) true-colour image of M83 shows strong 
unresolved emission both in the starburst nuclear region 
and along the spiral arms. Bottom panel: full-resolution 
true-colour image of the central region (size: 3' x 3'J. In 
both images, red, green and blue correspond to emission 
m the 0.3-1.0 keV, 1.0-2.0 keV and 2.0-8.0 keV bands 
respectively. Both images have been adaptively smoothed. 



3. Emission from the nuclear region 

The Chandra data reveal that M83 has a very highly struc- 
tured nuclear region (Figure 2, bottom panel). Fifteen dis- 
crete sources are detected within a radius of w 16" (« 
290 pc) from the centre of symmetry of the outer optical 
isophotes. A s pectral analysis of t he two brightest sources 
is presented in Soria fc Wu (2002) , and an analysis of other 
point sources will be presented in Soria et al. (in prepara- 



tion). We removed the point sources and extracted counts 
from concentric annuli to construct radial brightness pro- 
files of the unresolved emission in the nuclear region. We 
found that the brightness is approximately constant in 
a circular region up to a radius of 7" and then declines 
radially with a power-law like profile. The azimuthally- 
averaged profile is well fitted by a King profile, with a 
core radius of 6"7 ± 0"5 (« 120 pc), and a power law with 
a slope of —1.9 ± 0.2 beyond the core. 
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Table 1. Luminosity of emission from discrete and unre- 
solved sources in the nuclear region (0.3-8.0 keV band). 
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Figure 3. The cumulative luminosity distributions of 
sources found in the 0.3-8.0 keV band is different for the 
population inside the 60" inner circle (which includes the 
starburst nuclear region and the main bar) and for the 
population outside of it (red and blue curves respectively). 
We interpret the single power-law distribution of the nu- 
clear sources as evidence of continuous, ongoing star for- 
mation. The relative scarcity of bright sources in the disk 
population may be the effect of aging, which suggests that 
star formation is less active there. See Section 2 for the 
conversion from counts to luminosity. 



luminosity (xlO 38 erg s 1 ) 


inside 7" 


inside 16" 


discrete sources 


7.7 


12.3 


unresolved sources 


8.0 


11.5 


optically-thin thermal component 


5.5 


7.8 


power-law component 


2.5 


3.7 


total 


15.7 


23.8 



optically thin gas, and emission from unresolved point- 
like sources (eg, faint X-ray binaries). Assuming that the 
latter contribution is responsible for the power-law com- 
ponent in the spectrum of the unresolved emission, we 
estimate that emission from truly diffuse thermal plasma 
contributes w 35% of the total luminosity (Table 1). 

Extrapolating the log N(>S) - log S curve for the nu- 
clear sources gives us another way of estimating the rela- 
tive contribution to the unresolved emission of truly dif- 
fuse gas and faint point-like sources. We obtain that un- 
resolved point-like X-ray sources inside 16" would have a 

10 38 



. ,. r . it- total luminosity of ~ 2.7 x 10 38 erg s (see Soria & Wu 
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sion i nside the 



imiei 1U" ciicre, excluding btre resolved 

point sources, and we fitted it using an absorbed, single- 
temperature vmekal plus power-law model. Assuming so- 
lar abundances, we obtain a best-fit electron temperature 
kT = (0.60^003) keV and a power-law photon index V = 
3.1^q'2- The predicted lines are not strong enough to ac- 
count for the data, leading to poor fit statistics (xt = 1.42, 
114 dof). Increasing the abundance of all the metals by the 
same constant factor does not improve the fit. 

We then assumed a different set of abundances, higher 
than solar for C, Ne, Mg, Si an d S, and slightly under- 
ubandant for Fe (see Table 1 in |5oria fc Wu 2002|) . This 
is physically justified if the interstellar medium has been 
enriched by type-II supernova ejecta and winds from very 
massive, young stars. We obtain a best-fit temperature 
kT = {0.58±Q °l) keV, and power-law photon index T = 
2.7t u j (xl = 0.99, 114 dof). 

We also estimated the total (resolved plus unresolved) 
luminosity from the circular regions within radii of 7" (this 
is approximately the region inside the outer dust ring) 
and 16" from the geometric centre of the X-ray emis- 
sion, using an absorbed, optically-thin thermal plasma 
plus power-law model. The total emitted luminosity in the 
0.3-8.0 keV band is « 15.7 x 10 38 erg s" 1 inside 7" and 
« 23.8 x 10 38 erg s _1 inside 16" (Table 1). Discrete sources 
contribute « 50% of the total luminosity. The unresolved 
emission is itself the sum of truly diffuse emission from 



2002 for details). Another possible contribution to the un 
resolved emission comes from photons emitted by the re- 
solved sources but falling outside the extraction regions, in 
the wings of the PSF. We estimate that this contribution 
is ^ 1.5 x 10 38 erg s _1 . Thus, the combined contribution of 
faint X-ray sources and emission in the wings of the PSF 
can account for the luminosity of the power-law compo- 
nent inferred from the spectral fitting of the unresolved 
emission. This also confirms that a substantial proportion 
(« 70%) of the unresolved emission is indeed due to truly 
diffuse gas rather than faint point-like sources. 

Diffuse X-ray emission is also clearly observed along 
the spiral arms (Figure 2, top panel). We compared the 
spectrum of the unresolved emission in the nuclear region 
and along the arms, and found that the unresolved arm 
emission has a thermal component dominating at lower en- 
ergies, and a power-law-like component, dominating above 
1 keV. The thermal component in the arms has a temper- 
ature T = 0.3lt^3 keV, much cooler than that in the 
nucleus. The power-law component has a photon index 
r = 1.42+q'q 8 . A more detailed analysis of the unresolved 
arm emission will be presented in Soria et al. (in prepara- 
tion) . 

4. Discussion 

Separating the discrete sources inside and outside a 60" cen- 
tral region reveals that the two populations have different 
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Figure 4- The spectrum of the unresolved emission in 
the nuclear region (inside a 16" circle) shows very 
strong emission lines from MgXI (1.33-1.35 keV), SiXIII 
(1.84-1.87 keV) and SXV (2.43 keV). The other ma- 
jor contributions to the spectrum come from C VI (0.37- 
0.44 keV), OVII (triplet, 0.56-0.57 keV), NelX (triplet, 
0.91-0.92 keV), NeX (1.02 keV), and the FeXVII L-line 
complex (0.6-0.9 keV). 



cumulative luminosity distributions. The log N(>S) - log 
S curve of the sources outside this radius (i.e., the disk 
population) shows a kink at a luminosity ks 6 x 10 37 erg s _1 
in the 0.3-8.0 keV band. No kink is seen for the sources 
inside the 60" radius (i.e., those located in the nuclear re- 
gion and along the bar). The slope of the log N(>S) - log 
S curve at its high-luminosity end is flatter for the nuclear 
population, implying a larger proportion of bright sources 
(possible black-hole candidates). We interpret this as ev- 
idence of a past star formation episode in the disk, while 
there is continuo us, ongoing star formation in the nuclear 
region ( |Wu 200l| ). 

About 50% of the total emission in the nuclear re- 
gion belongs to resolved discrete sources. We estimate that 
« 70% of the unresolved emission (35% of the total) is due 
to truly diffuse plasma, with the rest (15% of the total) 
coming from faint, unresolved point-like sources and pho- 
tons in the wings of the PSF outside the detection cells of 
the resolved sources. 



The X-ray spectrum of the unresolved nuclear compo- 
nent shows strong emission lines, and can be modelled as 
emission from optically-thin thermal plasma at kT « 0.6 
keV. Above-solar abundances of Ne, Mg, Si and S are re- 
quired to fit the spectrum, while Fe appears to be un- 
derabundant. This suggests that the interstellar medium 
in the starburst nuclear region has been enriched by the 
ejecta of type-II supernova explosions. Moreover, a high 
abundance of C and a high C/O abundance ratio can 
be the effect of radiatively-driven winds from metal-rich 
Wolf-Rayet stars with M > 40 M dGustafsson et al. 1999| ). 
Both effects are likely to be present in the nuclear region. 




channel energy (keV) 

Figure 5. The ratio between the unresolved emission in the 
spiral arms and the best-fit model to the unresolved emis- 
sion in the nuclear region shows that: a) the thermal com- 
ponent in the spiral arms is cooler (T « 0.31 keV in the 
arms, T 0.58 keV in the nucleus); b) a power-law com- 
ponent dominates the spiral-arm emission at energies ft 1 
keV. The two spectra have been normalised to the same 
number of counts. 



Strong unresolved emission is also detected along the 
arms. It is well fitted by a thermal component at kT ss 0.3 
keV and a power-law component (r « 1.5) dominating at 
energies 2j 1 keV. A study of this higher-energy component 
and a comparison with the Galactic Ridge emission will 
be presented in a work now in preparation. 
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